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The effect of antiferromagnetic spin fluctuations on two-dimensional quarter-filled systems is stud- 
ied theoretically. An effective t — J' — V model on a square lattice which accounts for checkerboard 
charge fluctuations and next-nearest-neighbors antiferromagnetic spin fluctuations is considered. 
From calculations based on large- N theory on this model it is found that the exchange interaction, 
J', increases the attraction between electrons in the d xy channel only, so that both charge and spin 
fluctuations work cooperatively to produce d xy pairing. 
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I. INTRODUCTION 

The BEDT-TTF (bis-ethyleneditiotetrathiafulvalene) 
family of quarter-filled layered organic materials with 
the 6, a and j3" arrangements of the molecules dis- 
play a subtle competition between metallic, charge or- 
dered insulating and superconducting phases 1 ^. They 
are examples of strongly correlated electron systems for 
which their electronic states are theoretically described 
by a 2D extended Hubbard model at 3 /4- filling of elec- 
trons (1/4-filling of holes) for the HOMO of the BEDT- 
TTF molecules^. The nearest-neighbors intermolccular 
Coulomb interaction, V , is a crucial ingredient as, at one- 
quarter filling, the on site Coulomb repulsion, U, by itself 
cannot describe charge ordering phenomena 4 . The ex- 
tended Hubbard model at this filling has been previously 
studied— through large- N and slave-boson approaches in 
the [/-infinite limit as well as with exact diagonalization 
on small clusters^ at finite-t/. Several issues related to 
charge ordering phenomena have been addressed. A tran- 
sition from a metal to a checkerboard charge ordered in- 
sulating state at a finite V — V c has been found. Close to 
this charge ordered phase, superconductivity in the d xy 
channel appears- induced by strong charge fluctuations. 
Dynamical 7 properties of the metallic phase in the pres- 
ence of short range charge fluctuations have been found 
to be anomalous in agreement with experimental data 8 . 

Large- N methods and slave bosons are useful for the 
study of the effect of charge fluctuations on various elec- 
tronic properties as they can be included at 0(1/N), 
however, spin fluctuations are typically neglected un- 
less complicated 0(1/ N 2 ) contributions are considered. 
The effect of spin fluctuations on superconductivity in 
quarter-filled systems has been addressed through RPA 
(Random Phase Approximation) calculations which con- 
sider antiferromagnetic instabilities induced by the on- 
site Coulomb interaction, U 9], finding that d xy super- 
conductivity still prevails. However, well in the charge 
ordered insulating phase it is known that the spins or- 



der antifcrromagnctically due to the presence of a next- 
nearest neighbor spin exchange interaction^. This spin 
interaction results from a 'ring' exchange process appear- 
ing at fourth order in t and acts along the diagonals of 
the square lattice reading: J' = 4t 4 /9V 3 in the U — > oo, 
and V » t limits. Exact diagonalization on 16-site clus- 
ters indicate that the (ir/2,ir/2) antiferromagnetic spin 
arrangement follows closely the (tt, tt) checkerboard ar- 
rangement of the charge^ as V/t is increased from the 
metal to the charge ordered phase. These results suggest 
that remnants of the exchange interaction, J', generated 
in the insulating phase can survive in the metallic phase 
where short range charge ordering is present. It is then 
the purpose of the present work to analyze the influence 
of this exchange coupling, J', on the superconducting 
instabilities previously found 6 induced by charge fluctu- 
ations. As this J' acts along the diagonals of the lattice it 
is conceivable that similarly to the d x 2_ y 2 superconduc- 
tivity appearing in the t — J model close to half-filling 
induced by antiferromagnetic spin fluctuations, the J' 
appearing in the quarter-filled t — J' — V model can in- 
duce d xy pairing. 

The paper is organized as follows. In Sect. [H]we in- 
troduce the quarter-filled t — J' — V model and provide a 
phase diagram obtained from the large- N approach used 
here. In Sect. IIIII we discuss the superconducting phase 
and study in detail the effect of J' on the pairing sym- 
metry. Finally in Sect. II VI we summarize our results and 
point out their connection to the electronic properties of 
quarter-filled layered organic molecular crystals. 



II. THE t-J'-V MODEL AND ITS PHASE 
DIAGRAM 

In order to explore the above possibility an effective 
quarter-filled t — J' — V model is introduced as a natu- 
ral extension of the extended Hubbard model previously 
studied (see Fig. ^for a schematic sketch of the effective 
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interactions). The model reads 
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FIG. 1: A square lattice with a nearest-neighbors Coulomb re- 
pulsion V and hopping amplitude t and a next-nearest neigh- 
bors spin exchange coupling J'. 



where £y and Vij are the nearest-neighbors hopping and 
Coulomb repulsion parameters, respectively, between 
sites i and j on the square lattice. J?, is the antiferro- 
magnetic exchange interaction between second-neighbors 
sites, cf ia and Ci„ are the fermionic creation and destruc- 
tion operators respectively under the constraint that dou- 
ble occupancies of lattice sites i are excluded. Si and 
are the spin and the fermionic density, respectively. 

The model Q is studied by using the large- N approach 
for Hubbard operators^ recently extended to the case of 
finite (see also Ref.Q for the J = case). In the ap- 
pendix we give details about this approach. The method 
has been thoroughly tested on the t — U — V model by 
comparing dynamical properties with exact diagonaliza- 
tion calculations on small clusters^, finding good agree- 
ment for the behavior of charge and spectral functions 
as well as spectral densities close to the charge ordering 
transition. This agreement can be attributed to the fact 
that the infinite- U limit has been considered and neglect- 
ing the nearest-neighbor exchange J is justified since, at 
finite V, the charge tries to sit in every other site of the 
square lattice making J ineffective. 

Although the t— J'— V^-model is only justified for values 
of V sufficiently large, we have explored the full param- 
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FIG. 2: Phase diagram obtained from the large-iV approach 
in the V — J' plane. The solid line represents V c i. e. the crit- 
ical line signalling the onset of checkerboard charge ordering 
(CO). Between the solid and dashed lines d xy superconduc- 
tivity (SC) appears. The region below the dashed line corre- 
sponds to the metallic (M) phase. Superconductivity is found 
to be more robust closer to the CO line and for larger J'. 



eter range for completeness. In the t — J' — V model, 
the J' is dynamically generated when the charge is or- 
dered within the checkerboard pattern through a 'ring' 
exchange process. Hence, the J' becomes effective only 
when some sort of checkerboard charge ordering is al- 
ready present in the system (either short or long range 
charge order). This means that the system should be 
charge ordered or sufficiently close to the charge ordering 
transition for the t — J' ~ V model to be meaningful. The 
situation is different in the t — J model for the high-T c as, 
in this case, the J is generated through a super-exchange 
process at large on-site U. 

A full phase diagram summarizing our results obtained 
from large- N theory on the t — J' — V model is shown in 
Fig. |21 where metallic (M), charge ordered (CO) and su- 
perconducting (SC) phases occur. We start with a discus- 
sion of the charge ordering transition. The critical value, 
V c , signalling the charge ordering (CO) transition of the 
metallic phase as obtained from the divergence of the 
static charge susceptibility (see appendix), is displayed 
as a solid line in Fig. [2J For J' = 0, the system charge 
orders at V = V c ~ 0.65t, as previously found^. The 
value of V c is found to decrease with increasing J' which 
can be easily understood from the following. By increas- 
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FIG. 3: Behavior of the effective potential between quasipar- 
ticles, V eff {q, q), for J' = 0, J' = O.li and J' = 0.2t, for fixed 
= 0.5t. For this parameters the system is within the metal- 
lic phase close to the charge ordering transition (see phase di- 
agram in Fig. As J' is increased, V e ff(q, q) becomes more 
anisotropic and more attractive near the momentum transfer 
q ~ (7r,7r). Divergences appearing at (n,n) for increasing J' 
favor superconductivity in the d xy channel. 



ing J', checkerboard charge ordering is favored because 
two electrons are antiferromagnetically attracted when 
they are sitting along the diagonals of the square lattice 
so that a smaller V is effectively needed to induce CO. 
In a similar way, if we switch on J' for a fixed V < V c , 
the CO phase is favored. The dashed-dotted line in Fig. 
121 marks the onset for superconductivity. Between solid 
and dashed lines the superconducting effective coupling 
in the d xy channel (see below) becomes negative indica- 
tive of superconducting pairing. Superconductivity in 
the d xy channel is found to be more robust for larger J' 
and closer to CO, as one would expect as in either way a 
stronger attraction is felt between quasiparticles in every 
other site of the lattice. 



III. SUPERCONDUCTIVITY IN THE t-J'-V 
MODEL 

In order to understand the model proposed we first dis- 
cuss related large- N studies performed in the well known 
models such as the t — J model relevant to the cuprates. 
The case of a nearest-neighbors exchange coupling, J, 
has been studied using a Slave-boson techniques have 
been used in combination with al/N expansion to study 



superconducting instabilities on t — J modeU^ and Hub- 
bard modeiii Both the present approach (as described 
in the appendix) and the slave-boson one, to O(l), lead 
to fermions renormalized by the presence of the Coulomb 
interaction. Indeed, superconductivity in the model can 
only appear at 0(1/N). By adding an extra J' term 
to the bare extended Hubbard hamiltonian (t — J' — V- 
model) we can treat both J' and F-terms at the same 
level of approximation, i. e. through 0(1/N) and ana- 
lyze possible superconducting instabilities. A closely re- 
lated work by Vojta— analyzes recently coexistence of su- 
perconductivity and checkerboard charge ordering within 
a t — J — V model close to half-filling relevant to STM 
experiments on the cuprates^. 

Superconductivity is then investigated, within our 
model, by calculating the effective interaction, V4//(q), 
through 0(1/N) (see appendix), between fermions for 
finite J' at one-quarter filling. The potential, 14// (q) 
is plotted along the q = (q, q) direction in Fig. [31 as 
a function of q for different J': J' = 0, J' = Q.lt and 
J' = 0.2t, and fixed V — 0.5t. For this set of parameters 
the system is always metallic but close to the charge or- 
dering transition as can be seen from Fig. [21 For J' = 0, 
V4//(q) is repulsive and anisotropic as found previously 
(see Fig. 13 of Ref.0). As J' is increased, V e ff(q,q) 
becomes more anisotropic and more attractive near the 
momentum transfer q ~ (jr, n). This behavior favors su- 
perconductivity in the d xy channel as the J' attracts the 
charge tending to form the checkerboard pattern along 
the diagonals of the square lattice. 

There are two kinds of interactions contributing to 
V e ff (q). For J' = 0, the effective interaction close to CO 
is mainly dominated by charge fluctuations associated 
with the collective excitations near (n, 7r)l. Therefore, 
this kind of pairing interaction, like the phonon mecha- 
nism in simple metals, is mainly retarded and occurs in 
momentum space. For V — 0, the main effective interac- 
tion is of magnetic origin, unretarded and short range in 
real space. For J' and V different from zero both kind 
of interactions contribute cooperatively to the binding 
energy of the Cooper pairs. 

We use this effective potential to compute the effec- 
tive couplings in the different pairing channels or irre- 
ducible representations of the order parameter, i (i = 
(d x 2_ y 2, d xy ,p)). In this way we can project out the in- 
teraction with a certain symmetry. The critical temper- 
atures, T c , can then be estimated in weak coupling from: 
T C i = 1.13a;o exp(— l/|Aj|), where loq is a suitable cutoff 
frequency and Xi are the effective couplings with different 
symmetries. These are defined aai 
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where the functions <?i(k), encode the different pairing symmetries, and i>k are the quasiparticle velocities at 
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FIG. 4: The dimensionless superconducting coupling A in 
the p (X p ) (dashed line), d x 2_ y 2 (Xd 2^2) (dashed-dotted line) 
and d xy (Ad XH ) (solid line) channels as a function of J' for a 
given value of V = 0.5t in the metallic phase close to the 
CO instability. The second-neighbors antiferromagnetic ex- 
change, J', favors superconductivity in the d xy channel only. 

the Fermi surface. The integrations are restricted to the 
Fermi surface. A^ measures the strength of the interac- 
tion between electrons at the Fermi surface in a given 
symmetry channel i. If A,; > 0, electrons are repelled. 
Hence, superconductivity is only possible when Ai < 0. 

The coupling strength A for d xy superconductivity has 
been found£*£ to be very small (see Fig. 15 of Ref.Q) 
previously. For these small couplings, the corresponding 
superconducting critical temperature T c is predicted to 
be very low. Other indications for pairing come from 
exact diagonalization calculations of the binding energy 
of two holes which becomes negative near the charge or- 
dering transition*. In Fig. 0] we present results for the 
dimensionless superconducting coupling A in the p (A p ), 
d x 2_ y 2 (Ad x 2_ 2 ) an d d xy {Xd xy ) channels as a function of 
J' for a given value of V = Q.5t close to the CO instabil- 
ity. 

Clearly, superconductivity becomes more favorable in 
the d xy channel only as J' is increased. The influence of 
J' on the da;2_j,2-symmetry (dashed-dotted line in Fig. 

is rather weak. When J' increases, superconducting 
couplings are more repulsive in the p-channel and even 
more in the s-symmetry (not shown) channel. 

Note also the large difference between the A values for 
J' = C 7 and J' = 0.3t. For J' = 0.3i, A is one order of 
magnitude larger than for J' = at V = 0.5i. This be- 
havior, which only occurs in the d xy channel, shows the 
strong influence of the second neighbors effective antifer- 
romagnetic exchange coupling J' on d xy superconductiv- 
ity. 

Although we find a substantial enhancement of pairing 
in the d xy channel with J', the associated T c is yet very 
low from simple estimates: T c ~ e _1 /' Ai ' ~ 6 ■ 10~ 8 t, 
which is tiny even taking the most favorable case of J' = 
0.3i (see Fig. 0} for pairing, for which, Xd ly ~ —0.06. 



Small T c values have been also found by Motrunich and 
Lee-£ in the context of NazCo02, although their cou- 
plings are typically larger than ours because, in the case 
of Na^Co02, several sections of the Fermi surface are con- 
nected by the charge ordering wavevector for x = 1/3. It 
is worth noting that in our approach we have not consid- 
ered the renormalization of the quasiparticles, which en- 
hances the effective mass, that occurs close to the charge 
ordering transition due to V . This could be taken into 
account by including self-energy effects in the calculation. 
One would then have to consider that the hopping am- 
plitudes are effectively renormalized by the quasiparticlc 
weight, Z, which decreases near the charge instability 7 . 
This would transform the bare hopping, t, to t e ff = Zt. 
Hence, the ratios J/t e ff and V/t e ff would be enhanced, 
effectively increasing |Ad x | which would, in turn, signif- 
icantly increase the estimated T c . 



IV. CONCLUSIONS 

In conclusion, we have found that superconductivity 
with d xy symmetry induced by charge fluctuations is 
strengthened by antiferromagnetic spin fluctuations in- 
duced by the exchange coupling, J', in a t — J' — V 
model for quarter- filled systems. This can be intuitively 
understood by analogy with the more standard case of 
d x 2_ y 2 superconductivity induced by the superexchangc 
coupling J in the t — J model if one realizes that the J' 
acts in directions rotated by 45° with respect to J. Both 
charge and spin fluctuations are then found to work to- 
gether cooperatively to produce d xy superconductivity in 
the quarter-filled t — J' — V model. 

The way superconductivity behaves in the t — J' — V 
model proposed here, could be viewed as a two step pro- 
cess in which the charge fluctuations are responsible for 
the onset of SC in the first place and subsequently the 
dynamically generated spin exchange coupling J' would 
strengthen the binding between electrons forming the 
Cooper pairs. 

Critical temperatures are found to be too small com- 
pared to experimental values (which are of the order of a 
few Kelvin). The smallness of T c is related to the small 
Fermi surface associated with the one-quarter filling of 
the system. Due to this fact, there are no two points at 
the Fermi surface connected by the (ir, if) CO wavevec- 
tor which makes the interaction less effective in produc- 
ing Cooper pairs than in nearly antiferromagnetic metals 
close to half-filling. A similar situation arises in the t — J 
model— at one-quarter filling in which superconductivity 
in the d x 2_ y 2 channel is found although with rather small 
attractive effective couplings. Only at dopings close to 
half-filling (for doping levels of at most 0.15-0.2) the cou- 
plings are found to be substantial— This is because, 
in this case, larger Fermi surface sections are effectively 
connected by the AF (tt, it) wavevector. 

An important finding that derives from our work is 
that including J' is crucial in order to enhance T c suf- 
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ficiently. If only the V is taken into account, the values 
of obtained would be astronomically small as pre- 
viously notedi. A more sophisticated theory including 
the renormalization of the quasiparticles could enhance 
the estimates for T c even further. Measured effective 
masses in quarter-filled layered organics correspond to 
values of about: m* /m =1 — 2, leading to J' jt larger by 
a factor of 2. Considering the dependence of A<j with 
J' ft shown in Fig. 0] this many-body effect can lead to 
large enhancements in T c as T c depends exponentially 
with Ad • Hence, T c 's of a fraction of a Kelvin can be 
obtained in the most favorable case, for V — O.bt and 
J/t = 0.3*, considering t ~ 0.1 e\£. 

It has been recently suggested by Coldea et. al. 
Ref. 19] that certain experimental observations are 
consistent with the charge mediated superconductiv- 
ity scenario. For instance, the unit cell volume 
of /3"-(BEDT-TTF) 4 [(H30)M(C 2 1 )3]- Y [H and a- 
(BEDT-TTF) 2 MHg(SCN) 4 salts |9, is found to in- 
crease from metallic to insulating salts by changing 
M and Y. Superconducting salts such as a-(BEDT- 
TTF) 2 NH 4 Hg(SCN) 4 and the recently analyzed 0"- 
(BEDT-TTF) 4 [(H30)Ga(C204)3]-C6H 5 N0 2 are found 
to have unit cell volumes right between their respective 
metallic and insulating salts. Increasing the unit-cell vol- 
ume is translated to an increase in V/t as well as in J'/t 
within our model which drives the system closer to the 
charge ordered state. In the critical region between the 
metal and charge ordered phase, V < V c , superconduc- 
tivity is predicted to appear—. 

However, a definitive test for unconventional d-wave 
pairing in quasi-two-dimensional quarter-filled organics 
is yet missing. Possible experimental probes could come 



The five Hubbard ^-operators X? a ' and X°° are 
boson-like and the four Hubbard A-operators Xf and 
X® a are fermion-like. The names fermion-like and boson- 
like come from the fact that Hubbard operators do 
not verify the usual fermionic and bosonic commutation 
relations^. 

From the above relations we note that Xf° — c] ia and 

The Hubbard operators satisfy: 
a)the completeness condition 

X°° + Y,Xr = h (A2) 



from measurements of the Knight shift and NMR relax- 
ation rate in the metallic phase close to the CO transi- 
tion. In contrast to the k-(BEDT-TTF) 2 X superconduc- 
tors, in which large enhancements of the Korringa ratio 
are found due to their closeness to a Mott phase, in the 
quarter-filled systems studied here there should be no 
enhancement of the Korringa ratio. Finally, the depen- 
dence of T c on impurities and disorder can also be used to 
distinguish d-wave superconductivity from conventional 
s-wave pairing as recently pointed oufcSi. 
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APPENDIX A: APPENDIX: LARGE- iV 
APPROACH FOR t-J'-V MODEL 

In this appendix, we will give a summary of the path 
integral large-N approach for Hubbard operatora£ii*i£ 
used in the present paper. More detailed version of the 
formalism including the exchange coupling can be found 
ini£. 

First we introduce Hubbard operators which are re- 
lated with the usual fermionic operators by 



(Al) 



I 

which is equivalent to imposing that " double occupancy 
at each site is forbidden. 

b)the commutation rules 



[Xf, X] 5 ]± = 5ij {S^Xf ± 5 aS xf) (A3) 

where the + sign must be used when both operators are 
fermion-like, otherwise it corresponds the — sign. 

On the basis of Hubbard X-operators, the t — J' — V 
Hamiltonian (1) is of the form: 



J 



l»^(l-i ClJ )c„, Xf = {X°°)\ Xf =c\ aCia ,. 
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h(x)= {ti j xfx^+h.c.) + \ Ji 3 (xfx™ -x™xf)+ Vijxrxf*' -»J2 x rm 



<i]>,rr 



<ij>;<7 



Our starting point is the path integral partition function Z written in the Euclidean form 



z = J vxf s[xf° + J2*r i] S[X? 



crO vOcr' 



xrx\ 

X?" 



-] x (sdetM^ s )f exp (- / dr L E (X,X)) 



(A5) 
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The Euclidean Lagrangian Le(X, X) in (A5) is 



L E {X,X) 



(xr x 



ctO 



X°° Xf a ) 



2 ^ 



x. 



00 



+ H(X$A6) 



In this path integral we associate Grassmann and 
usual bosonic variables with Fermi-like and boson-like 
A-operators, respectively. 

It is worth noting at this point that the path integral 
representation of the partition function (A5), looks dif- 
ferent to that usually found in other solid state problems. 

The measure of the integral contains additional con- 

i 

straints as well as a determinant, (sdetMAB)t ■ Also the 
kinetic term of the Lagrangian (A6) is non-polynomial. 
The determinant reads 



(sdetM AB )* = 1/ 



(-A 00 ) 2 ' 



(A7) 



and is formed by all the constraints of the theory 
We now discuss the main steps needed to introduce a 
large- N expansion of the partition function (A5). First, 
we integrate over the boson variables A CTCT using the sec- 
ond S- function in (A5). We extend the spin index a = ±, 
to a new index p running from 1 to A. In order to get a 
finite theory in the N — ► oo limit, we re-scale tij to tij/N, 
Vij to Vij/N and J[- to J-j/N. The completeness condi- 
tion is enforced by exponentiating Xf + ^ Xf p = N/2, 



with the help of Lagrangian multipliers Aj . We write the 
boson fields in terms of static mean-field values, (r , A ) 
and dynamic fluctuations 



= Aro(l + <5i?,) 

A, = X a + SX t , (A8) 
and, we make the following change of variables 



I ip 



1 



i 



X 



pO 



(A9) 



where and fi p are Grassmann variables. 

The exchange interactions can be decoupled in terms of 
the bond variable A^ through a Hubbard-Stratonovich 
transformation, where A^- is the field associated with the 



We write the An fields in 



quantity V — . 

term of static mean field values and dynamics fluctua- 
tions A^ = A(l + r) t + iAf), where 77 can take two val- 
ues associated with the bond directions 771 = (1,1) and 
rj2 = (—1,1) in real space. 

Introducing the above change of variables into Eq. 
(A5) and, after expanding the denominators 1/(1 + 5R), 
we arrive at the following effective Lagrangian: 
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Feynman Rules 

Propagators 




FIG. 5: Summary of the Feynman rules. Solid line represents 
the propagator C?( ) f° r the correlated fermion X 0a . Dashed 
line represents the 6x6 boson propagator D( ) for the 6- 
component field 8X a . The component (1, 1) of this propaga- 
tor is directly associated with the X 00 charge operator. Doted 
line is the propagator B for the boson ghost field Z p . A pp and 
^ab represent the interaction between two fermions f p and 

one and two bosons 5X a respectively. F pp and r£{| represent 
the interaction between two ghost fields Z v and one and two 
bosons 8X a respectively. 



Leff = - 2 E {M + ftpfip) (1 - SRi + SRf) + ]T farof&fip + h - c -) ~ M E ~ 6R * + SR ^ 

i 7 p <ij>,p i,p 

^ O AT 

+ Nr a J2 6X > 5R i + E ftjiv^ SRi) SX t + —A 2 £ [(rf) 2 + (A?f] 

i i,p if] 

- A £ (/+/ jy +/+j ip ) [1-1(^ + ^0 + ^^ + 1(^ + ^)1 

<ij>,P,p' 

<ij>;P;P' 

+ Nrl (Vij ~ \j i j)5R i 8R j - £ Z+, (l - «fc + 5i? 2 ) Z ip , (A10) 



where we have changed to fi — Xq and dropped constant 
and linear terms in the fields. 

The last term of (A10) results from the path inte- 
gral representation of the determinant which uses N- 
component boson ghost field Zjjiii*^. 

Looking at the effective Lagrangian (A10), the Feyn- 
man rules can be obtained as usual. The bilinear parts 
give rise to the propagators and the remaining pieces are 
represented by vertices. Besides, we assume the equation 
(A10) written in the momentum space once the Fourier 
transformation was performed. 

To leading order of l/N, we associate with the N- 



r 



component fermion field f p , connecting two generic com- 
ponents p and p', the propagator 



G(o)p P '(k, v n ) 



Opp' 



iv n - (E k - fi) 



(All) 



which is O(l) and where E k = —2tro(cosk x + cosk y ) — 
2Acosk x cosk y , is the electronic dispersion to leading or- 
der. 

The quantities k and v n are the momentum and the 
fermionic Matsubara frequency of the fermionic field, re- 
spectively. 
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a) Irreducible boson self-energy 

n = n (1) + n (2) + n® + n <4) 

a^_b = £_Q_b + aQb + i-'>- b + tL/J> 

b) Effective interaction between fermions 

V eff = | 

FIG. 6: a) The four different contributions (i = 1, 2, 3, 4) 
to the irreducible boson self-energy fl a f,. b) Effective inter- 
action between fermions. Only two three- legs vertices con- 
tribute. 



The mean field values r$ and A must be determined 
minimizing the leading order theory. From (A8) and the 
completeness condition, ro is equal to 6/2, where 5 is the 
hole doping away from half-filling. 

On the other hand, minimizing with respect to A we 
obtain A = 4pjj- cosk x cos k y np(E\^ — fi), where rip 
is the Fermi function and N s is the number of sites in the 
Brillouin zone (BZ). 

For a given doping, 6, the chemical potential \x and 
A must be determined self-consistently from (1 — 5) = 

We associate with the six component SX a = 
(5R , SX, r 1 , r 2 , A 1 , A 2 ) the inverse of the propagator 
(which is 0(1/ N)), connecting two generic components 
a and b, 



D 



(O)afc 



(q.O =N 



(% 


ro 














ro 
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±A 2 

















' 








V o 











±A 2 



(A12) 



where 7 q = r 2 [4V(cos(q x )+cos(qy))-AJ'cos(q x ) cos(q y )] 
The quantities q and u n are the momentum and the 
Bose Matsubara frequency of the boson field, respec- 
tively. 

We associate with the TV-component ghost field, Z p , 
the propagator connecting two generic components p and 
P', 

Bppi = Spp' , (A13) 

which is 0(1) 

The non-quadratic terms in (A10) define three and four 
leg vertices which are 0(l).(The full expressions for the 
vertices associated with the large- N approach are given 
in Ref.[l^|). In Fig. [3]the Feynman rules associated with 
the large- N approach are summarized. 

The charge-charge correlation function is directly as- 
sociated with the element (1,1) of the boson propagator 
D(o) a b(q, w„) (the inverse of Eq. (A12)). As in Ref.0, 
up to 0(1/ N), D (c[,oj n ) is renormalized to _D(q, w n ) by 
an infinite series of diagrams of 0(1/ N) (Fig. and 
reads: 

D-^cbWn) = Do X (cbWn) -n(cbw„), (AH) 

where II(q, cu n ) is the boson self-energy for which explicit 
expressions are given in Ref. ITU 



The superconducting effective interaction between 
fermions, V e ff(q,u} n ), can be calculated using the Feyn- 
man rules of Fig. [5] Fig. |H|d) represents the diagram 
involved in the calculation of 14// (q, w n ). The analyti- 
cal expression for this diagram is 14// = A a D a b^-b where 
D a b is the propagator of the bosonic field which contains 
the irreducible self energies of Fig. 0a and A a is the three 
leg vertex of Fig. [SJ Looking at the order of the propa- 
gators and vertices we see that V e //(q, w„) is 0(1/ N). 

To conclude this appendix we make contact with 
closely related approaches such as slave boson formula- 
tions. In contrast to slave boson theories: (a) Greens 
functions are calculated in terms of the original Hubbard 
operators, (b) fermions, fi P , appearing in the theory are 
proportional to the Fermi- like A-operator X op (see (A9)) 
to all orders in the 1/N expansion; not only to leading 
order—, (c) as our path integral is written in terms of 
A-operators we do not need to introduce a priori any 
decoupling scheme, and (d) ro is the mean value of A 00 
which is a real field associated with the number of holes 
(see Eq. (A8)) and not with the number of holons. At 
leading order (N — > oo or 0(1)) and V = 0, our formal- 
ism is equivalent to slave boson approaches. However, 
at the next to leading order (0(1/N)), (which is neces- 
sary to calculate one-electron properties such as the elec- 
tron self energy E(k, ui) and the electron spectral func- 
tion A(k, w)), the two formulations do not coincide. The 
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differences between the two formulations are not yet com- 
pletely established. Our theory has the advantage that it 
does not require the introduction of gauge fields like in 
slave boson approaches. Hence, through order 0(1/N) 
we do not need to take care of gauge fluctuations nor 
Bose condensation (note that Eq. (A8) does not mean 



Bose condensation). This is important because for the 
doped Hubbard model the gauge fluctuations are known 
to significantly change the physics. Careful numerical 
work will determine the improvements of the present ap- 
proach with respect to slave boson formulations. 
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